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Abstract This review discusses the rale of radiative heating and cooling, as well 
as self-gravity, in shaping the nature of the turbulence in the interstellar medium 
(ISM) of our galaxy. The ability of the gas to radiatively cool, while simultaneously 
being immersed in a radiative heat bath, causes it to be much more compressible 
■ than if it were adiabatic, and, in some regimes of density of density and tempera- 

ture, to become thermally unstable, and thus tend to spontaneously segregate into 
separate phases, one warm and diffuse, the other dense and cold. On the other hand, 
turbulence is an inherently mixing process, thus tending to replenish the density 
and temperature ranges that would be forbidden under thermal processes alone. 
The turbulence in the ionized ISM appears to be transonic (i.e, with Mach num- 
bers Ms ~ 1), and thus to behave essentially incompressibly. However, in the neutral 
medium, thermal instability causes the sound speed of the gas to fluctuate by up 
to factors of ^ 30, and thus the flow can be highly supersonic with respect to the 
dense, cold gas, although numerical simulations suggest that the supersonic veloc- 
ity dispersion corresponds more to the ensemble of cold clumps than to the clumps' 
internal velocity dispersion. Finally, coherent large-scale compressions in the warm 
neutral medium (induced by, say, the passage of spiral arms or by supernova shock 
waves) can produce large, dense clouds that are affected by their own self-gravity, 
and begin to contract gravitationally. Because they are populated by large-amplitude 
density fluctuations, whose local free-fall times can be significantly smaller than that 
of the whole cloud, the fluctuations terminate their collapse earlier, giving rise to a 
regime of hierarchical gravitational fragmentation, with small-scale collapses oc- 
curring within larger-scale ones. Thus, the "turbulence" in the cold, dense clouds 
may be dominated by a gravitationally contracting component at all scales. 
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1 Introduction 



Our galaxy, the Milky Way (or simply, the Galaxy) is a flattened conglomerate of 
stars, gas, dust, and other debris, such as planets, meteorites, etc., with a total mass ~ 
IO^^Mq, where Me = 2 x 10^^ g is the mass of the Sun. Most of this mass is believed 
to be in a roughly spherical dark matter halo, while ^ IO^'Mq are in stars, and 
are co ntained in the gaseous component, mostly confined to the Galactic 
disk dcoxllioool) . 

The gasesous component may be in either ionized, neutral atomic or neutral 
molecular forms, spanning a huge range of densities and temperatures, from the 
so-called hot ionized medium (HIM), with densities n ^ 10^^ cm^^ and temper- 
atures T ^ 10^ K, through the warm ionized and neutral (atomic) media (WIM 
and WNM, respectively), with n ^ 0.3 cm^^ and T ^ 10^ K and the cold neutral 
(atomic) medium (CNM, n ^ 30 cm^^, T ^ 100 K, to the giant molecular clouds 
(GMCs, n ^ 100 cm^^ and T ^ 10-20 K. GMCs can span several tens of parsecs 
across, and, in turn, have plenty of substructure, which are commonly referred to 



, size scales L of a few parsecs), clumps (n 



10* cm- 



as clouds (n ^ 10 cm 

L ^ 1 pc), and cores (n ^ 10^ cm^^, L 0.1 pc), alth ough the temp erature of all 
molecular gas is remarkably uniform, ~ 10-30 K (e . g . . lFerrierell200 1 b . 

The gaseous component, along with the dust, a cosmic-ray background, and a 
magnetic field of mean intensity of a few jUG constitute what we know as the 
interstellar medium (ISM). Thi s medium is in most cases well described by the 
fluid approximation (IShul 1 19921 Ch. 1). Moreover, the ISM is most certainly tur- 
bulent, as typical Re ynolds numbers in it are ve ry large. For example, in the cold 
ISM, Re 10^-10^ (lElmegreen & Scalol |2004[ Sec. 4.1). This is mostly due to 
the very large spatial scales involved in interstellar flows. Because the ISM's tem- 
perature varies so much from one type of region to an other, so does the sound 



speed, and the flow is oft en super- or trans-sonic (e.g., iHeiles & Troland L 12003 



Elmegreen & Scalol 2004 , and references therein). This implies that the flow is sig 



nificantly compressible, inducing significant density fluctuations (Sec. 13. 2b . 

In addition to being turbulent, the ISM is subject to a number of additional phys- 
ical processes, such as gravitational forces exerted by the stellar and dark matter 
components as well as by the ISM itself (i.e., by its self- gravity), magnetic fields, 
cooling by radiative microscopic processes, and radiative heating due both to nearby 
stellar sources as well as to diffuse background radiative fields. All of this adds up 
to make the ISM an extremely complex and dynamical medium. 

Finally, the self-gravity of the gas causes local (i.e., spatially intermittent) grav- 
itational collapse events, in which a certain gas parcel within a molecular cloud 
goes out of equilibrium between its self-gravity and all other forces that oppose it, 
undergoing an implosion that leads the gas density to increase by many orders of 
magnitude, and whose end product is a star or a group ("cluster") of stars. 

In this review, we focus on the interaction between turbulence and the effects of 
radiative heating and cooling, which effectively enhance the compressibility of the 
flow, and have a direct effect on the star formation process. The plan of the paper is 
as follows: in Sec.|2]we first review the effects that the net heating and cooUng have 
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on the effective equation of state of the flow and, in the case of thermafly unstable 
flows, on its tendency to spontaneously segregate in distinct phases. Next, in Sec. 
[3] we discuss a few basic notions about turbulence and the turbulent production of 
density fluctuations in the compressible case, to then discuss, in Sec.|4] the interplay 
between turbulence and the heating and cooling. In Sec. |5] we discuss the likely 
nature of turbulence in the diffuse (warm and hot) parts of the ISM, while in Sec. 
|6] we do the same for the dense, cold atomic and molecular clouds. We conclude 
in Sec. |7] with a summary and some final remarks. Due to space limitations, we do 
not discuss magnetic fields, although we refer the intere sted reader to t he reviews 
bv lVazquez-Semadeni et al. (2000b.) „Cho et al.k2003.) . .Elmegreen & ScalQ(.20o4 . 
and lMcKee & Ostriker](l2007b . 



2 ISM Thermodynamics: Thermal InstabUity 



The ISM extends essentially over the entire disk of the Galaxy and, when consider- 
ing a certain subregion of it, such as a cloud or cloud complex, it is necessary to re- 
alize that any such subregion constitutes an open system, whose interactions with its 
environment need to be taken into account. A fundamental form of interaction with 
the surroundings, besides dynamical interactions, is through the exchange of heat. 
Indeed, the ISM is permeated by a radiation field, due to the combined shine of the 
stars in the disk. Moreover, a bath of relativistic charged particles, mostly protons, 
known as cosmic rays, also exists in the ISM. The se are believed to be acceler ated in 
strong shocks produced by supernova explosions ( IBlandford & Eichlen,ll987h . Both 
UV radiation and cosmic rays provide heating and ionization sou rces for the ISM 



at large (see, e.g.. lDalgarno & McCravill972tlWolfire et aUll995h . Finally, violent 



events, such as supernova explosions, can locally heat their surroundings to very 
high temperatures, causing local bubbles of hot, million-degree ionized gas. 

On the other hand, the ISM can cool by emission from ions, atoms and molecules, 

and by thermal emission fr om dust grains jPalgarno & McCravill972tlSutherland & Dopita , 



19931: IWolfire et all Il995h . The balance between these radiative heating and cool- 



ing radiative processes, along with the heat due to mechanical work and thermal 
conductivity, determine the thermodynamic properties of the ISM. These proper- 
ties will be the focus of the present section, rather than the detailed microphysical 
processes that mediate the radia tive transfer, about which there exist plenty of ex- 
cellent books and reviews (e.g., Dalgarno & McCravl Il972t : [Sutherland & Dopita , 
19931: lOsterbrock & FerlaiidibOOd) . 

Globally, and as a first approximation, the ISM is roughly isobaric, as illustrated 
in the left panel of Fig. [T] As can be seen there, most types of regions, either dilute or 
dense, lie within an order of magnitude from a thermal pressure P 3000 K cm^^Q 
The largest deviations from this pressure uniformity are found in HII regions, which 
are the ionized regions around massive stars due to the star's UV radiation, and 



^ It is customary in Astrophysics to express pressure in units of [K cm ■']. Strictly speaking, this 
corresponds to P/k, where k is the Boltzmann constant. 
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molecular clouds, which, as we shall see in Sec. |6] are probably pressurized by 
gravitational compression. 





Fig. 1 Left: Thermal pressure in various types of interstellar regions. The points labeled coronal 
correspond essentially to what is referred to as the HIM in the text; intercloud regions refer to 
the WIM and WNM; diffuse to CNM clouds, and dark, globule molecular to molecular gas. From 
[Myers ( 197^. Right: Temperature dependence of the c ooling function. The labels i ndicate values 
of the ionization fraction (per number) of the gas. From lDalgamo & McCra\j jl972h . 



The peculiar thermodynamic behavior of the ISM is due to the functional forms 
of the radiative heating and cooling functions acting on it, which depend on the 
density, temperature, and chemical composition of the gas. T he right panel of Fig. Q] 
shows the temperature dependence of the cooling function A ( Dalgarno & McCravT 
1972h . 

To understand the effect of the functional form of the radiative heating and cool- 
ing functions on the net behavior of the gas, let us wr ite the cons ervation equation 
for the internal energy per unit mass, e. We have (e.g.. lShulll992h 



dt 



(7- l)eV-u + r-«A, 



(1) 



where we have neglected thermal conduction and heating by magnetic reconnection. 
In eq. ([U, u is the velocity vector, / is the ratio of specific heats, F is the heating rate 
per unit mass, n ^ p / pLnin is the number density, with /i the mean partucle mass 
and mn the hydrogen atom's mass. Note that nA is the cooling rate per unit mass. 
The first term on the right hand side is the P dV work per unit mass. 

In the simplest possible case, that of a hydrostatic (u = 0) and steady {d /dt — 0) 
state, eq. ([T]) reduces to the condition of thermal equilibrium. 
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r{p,T)^nAip,T), (2) 

where we have denoted explicitly the dependence of F and A on the density and 
temperature, and neglected the dependence on chemical composition. 

Let us now assume we have a gas parcel in thermal and hydrostatic equilibrium, 
at a temperature somewhere in the range where the slope of A with respect to T 
is negative (i.e., {dA/dT)p < O), 10^-^ < T < 10*' K (see the right panel of Fig. 
[Til. If we then consider a small isochoric (i.e., at constant density) increase in T 
of this fluid parcel, we see that A decreases. Because we had started from thermal 
equiUbrium, a drop in A impUes that now F > nA, a condition that increases the 
parcel's temperature even further, causing a runaway to higher temperatures, until 
the fluid parcel exits the temperature range where {dA/dT)p < 0, at T ^ 10^ K. 
This behavior is known as thermal instability (TI), and the condition (dA/dT)p < 
is known as the isochoric criterion for TI (Field 1965; see also the review by 
Vazquez-Semadeni et al. 2003). 

Another, less stringent condition for the development of TI is the so-called iso- 
baric criterion. This can be most easily understood as follows. Note that eq. (|2]i 
provides a relation between the density and temperature of the medium. This can be 
inserted in the equation of state for the gas0 allowing the elimination of the tem- 
perature from it and writing a barotropic equation of the form Peq = ^eq(p), where 
Peq is the thermal pressure under conditions of thermal equilibrium, and is a func- 
tion of the density only. Figure |2] shows the resulti ng density dependen ce of fgq for 



the atomic medium under "standard" conditions (IWolfire et al.L 119951) . The region 
above the graph of Peq(n) corresponds to nA > F, and therefore to net cooling. 
Conversely, the area under the curve corresponds to net heating, nA < F. From this 
figure, we observe that the slope of the graph of Peq vs. p is negative for densities in 
the range 0.6 ^n <5 cm^^ (i.e., —0.2 < log[o(«/ cm"^) < 0.7). 

Let us now consider a fluid parcel in this density range, and apply to it a small, 
quasi-static compression (i.e., volume reduction) to it, increasing its density. This 
displaces the parcel to the region above the thermal equilibrium curve, where net 
cooling occurs. Because the compression is quasi-static, the parcel has time to cool 
under the effect of the net cooling, which causes its pressure to decrease as it at- 
tempts to return to the thermal equilibrium curve. But thus it is now at a lower pres- 
sure than its surroundings, which then continue to further compress the parcel, and 
runaway compression sets in, until the parcel exits the regime where dPeq/dyi < 0. 
The density range where this condition is met is called the unstable range. 

Thus, if the mean density of the medium is in the unstable range, this mode of TI 
tends to cause the me dium to spontaneo usly segregate into a cold, dense phase and 
a warm, diffuse one dField et al. . Il969h . That is indeed the case of the ISM in the 



Galactic midplane in the Solar neighborhood, a fact which led Field et al. ( 19691) to 



propose the so-called two-phase model of the ISM: assuming dynamical and thermal 
equilibrium, the atomic ISM would consist of small dense, cold clumps (the CNM), 
immersed in a warm, diffuse background (the WNM). The clumps are expected to 



^ Given the low densities of the ISM, it is a good approximation to use the ideal gas equation of 
state. 
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Fig. 2 Thermal-equilibrium pressure Peq as a function of number density for "standard" conditions 
of metallicity and b ackground UV radia tion for the atomic medium. The horizontal axis gives 
logio(/Vcm-')- From lWolfire et al.Ul99^ . 



be small because the fastest growing mode of the instability occurs at vanishingly 
small scales in the absence of thermal conductivity, or at scales ^0.1 pc for th e 



Audit &Hennebelle, 2005). 



estimated thermal conductivity of the ISM (see, e.g 
Note that the isochoric criterion for Tl implies that the isobaric one is satisfied, 
but not the other way around. Technical and math ematical deta ils, as well as other 
m odes of TI , can b e found in the original paper by 'Field' and in the reviews 

bv lMeerson ( 19961) and bylVazquez-Semadeni et al. (2003^^ 

It is important to note that, even if the medium is not thermally unstable, the 
balance between heating and cooling implies a certain functional dependence of 
P<:q(p), which is often approximated by a polytropi c law of the form Peq °= p^" (e.g 



Vazquez-Semadeni et al.l Il99d) . where it should be noted that 



Elmegreeni 11991 : 

7e, which we refer to as the effective polytropic exponent, is not the ratio of spe- 
cific heats for the gas in this case, but rather a free parameter that depends on the 
functional forms of A and F. The isobaric mode of TI corresponds to Ye < 0. 



3 Compressible Turbulence 
3.1 Equations 

In the previous section we have discussed thermal aspects of the ISM, whose main 
dynamical effect is the segregation of the medium into the cold and warm phases. 
Let us now discuss dynamics. As was mentioned in Sec. [T] the ISM is in general 
highly turbulent, and therefore it is necessary to understand the interplay between 
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turbulence and th e effects of the net cooling (nA — F), which affects the compress 



ibiUty of the gas (iVazquez-Semadeni et al.Lll996h . 



The dynamics of the ISM are governed by the fluid equ ations, which, n eglect ing 
magnetic fields, comprise eq. O and (e.g.. lLandau & Lifshitz,.1959.:.Shuill992h 



<9u 



dt 



u Vu 



-pVu, 
VP 
P 

4-nGp 



vrv-u) 



(3) 

(4) 
(5) 



where (p is the gravitational potential, and v is the kinematic viscosity. Equation (|3]l 
represents mass conservation, and is also known as the continuity equation. Equa- 
tion (in is the momentum conservation, or Navier-Stokes equation per unit mass, 
with an additional source term representing the gravitational force V^/p. In turn, 
the gravitational potential is given by Poisson's equation, eq. (|5]l. Equations ([T]), 
(O, (IDi, and ^ are to be solved simultaneously, given some initial and boundary 
conditions. 

A brief discussion of the various terms in eq. (HJi is in order The second term on 
the left is known as the advective term, and represents the transport of i-momentum 
by the j component of the velocity, where / and j represent any two components of 
the velocity. It is responsible for mixing. The pressure gradient term (first term on the 
right-hand side [RHS]) in general acts to counteract pressure, and therefore density, 
gradients across the flow. Finally, the term in the brackets on the RHS, the viscous 
term, being of a diffusive nature, tends to erase velocity gradients, thus tending to 
produce a uniform flow. 



3.2 Reynolds and Mach numbers. Compressibility 



Turbulence develops in a flow when the ratio of the advective term to the viscous 
term becomes very large. That is. 



^ [u ■ Vu] 



V (v^u ~ 



V(V-u) 
3 



L 



U_ 



n -I 



UL 



= /?e> 1, 



(6) 



where is the Reynolds number, U and L are characteristic velocity and length 
scales for the flow, and 0' denotes "order of magnitude". This condition implies that 
the mixing action of the advective term overwhelms the velocity-smoothing action 
of the viscous term. 

On the other hand, noting that the advective and pressure gradient terms con- 
tribute comparably to the production of density fluctuations, we can write 



8 



E. Vazquez-Semadeni 



■Vu) 



&{VP/p) 
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Ap 
P 



(7) 
(8) 



where Ms = U /c^ is the sonic Mach number, and we have made the approxima- 
tion that AP/Ap ^ Cj, where Cs is the sound speed. Equation ([8]l then implies 
that strong compressibility requires Ms ^ 1. Conversely, flows with Ms <C 1 be- 
have incompressibly, even if they are gaseous. Such is the case, for example, of the 
Earth's atmosphere. In the incompressible limit, p = est., and thus eq. (|3]l reduces 
to V • u = 0. 

Finally, a trivial, but often overlooked, fact is that, in order to produce a density 
enhancement in a certain region of the flow, the velocity at that point must have a 
negative divergence (i.e., a convergence), as indicated by the continuity equation, 
eq. (O. It is very frequent to encounter in the literature discussions of pre-existing 
density enhancements ("clumps") in hydrostatic equilibrium. But it should be kept 
in mind that these can only exist in multi-phase media, where a dilute, warm phase 
can have the same pressure as a denser, but colder, clump. But even in this case, 
\ht formation of that clump must have initially involved the convergence of the flow 
towards the cloud, and the hydrostatic situation is applicable in the limit of very 
long times after the formation of the clump, when the convergence of the flow has 
subsided. 



3.3 Production of Density Fluctuations 

According to the previous discussion, a turbulent flow in which the velocity fluctu- 
ations are supersonic will naturally develop strong density fluctuations. Note, how- 
ever, that the nature of turbulent density fluctuations in a single-phase medium (such 
as, for example, a regular isothermal or adiabatic flow) is very different from that 
of the cloudlets formed by TI (cf. Sec.|2l). In a single-phase turbulent medium, tur- 
bulent density fluctuations must be transient, as a higher density generally conveys 
a higher pressureJl and therefore the fluctuations must re-expand after the compres- 
sion that produced them has subsided. 

For astrophysical purposes it is important to determine the distribution of these 
fluctuations, as they may constitute, or at least provide the seeds for, what we nor- 
mally refer to as "clouds" in the ISM. Because of the transient nature of turbu- 
lent density fluctuations in single-phase media, however, this distribution refers to 
a time-stationary population of fluctuations, although the fluctuations themselves 
will appear and disappear over times short compared to the time over which the 
distribution is considered. 



^ An exception would be a so-called B urgers' flow, which is characterized by the absence of the 
pressure gradient term jBurgerj,[T974h . 
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The probability density distribution (PDF) of the density field in turbulent isother- 
mal flows was first investig ated numerically, finding t hat, in the isothermal case, it 
posesses a lognormal form dVazquez-SemadeM, 19941). A theory for the emergenc e 
of this functional form was later proposed bv lPassot & Vazquez-Semadenil ( ll998h . 
in which the production of density fluctuations was assumed to arise from a suc- 
cession of compressive or expansive waves, each one acting on the value of the 
density left by the previous one. Because the medium contains a unique distribution 
of (compressible) velocity fluctuations, and because the density jumps in isothermal 
flow depend only on Mach number but not in the local density, the density fluctu- 
ations belong all to a unique distribution as well, yet each one can be considered 
independent of the others if the global time scales c onsidered are much lo nger than 
the autocorrelation time of the velocity divergence ( Blaisdell et al.l [l993h . Finally, 
because the density jumps are multiplicative in the density (cf. eq. [8]l, then they 
are additive in = Inp. Under these conditions, the Central Limit Theorem can 
be invoked for the increments in s, implying that s will be normally distributed. In 
consequence, p will have a lognormal PDF. 



In addition, iPassot & Vazquez-Semadenil (Il998h also argued that the variance of 
the density fluctuations should s cale linearly with Ms, a fact that has been investi- 
gated further by various grou ps ( Padoan et al. , 1997 ; Passot & Vazquez-Semadeiil 



1998t iFederrath et al.L l2008l) . In particular, using numerical simulations of com- 



pressible turbulence driven by either solenoidal (or "vortical") or compressible (or 
"potential") forces, the latter authors proposed that the variance of s is given by 



(7,= [\n{l +bM^ 



(9) 



where ^7 is a constant whose value depends on the nature of the forcing, taking the 
extreme values of b— 1/3 for purely solenoidal forcing, and b ~ I for purely com- 
pr essible forcing. The lognormal densi ty PDF for the one-dimensional simulations 
of lPassot & Vazquez-Semadenil (119981) . with its dependence on Ms, is illustrated in 
the left pa nel of Fig. [3] 



Finally, Passot & Vazquez-Semadenil ( 1998 . see also Padoan & Nordlund 1999) 
also investigated the case where the flow behaves as a polytrope with arbitrary values 
of 7e, by noting that in this case the sound speed is not constant, but rather depends 
on the density as Cs <^ p(Tc-i)/2^ implying that the local Mach number of a fluid 
parcel now depends on the local density besides its dependence on the value of 
the flow velocity. I ntroducing this dependence of Ms o n p in the expression for 
the lognormal PDF, IPassot & Vazquez-Semadenil (Il998l) concluded that the density 
PDF should develop a power-law tail, at high densities when 7e < 1, and at low 
densities when 7e > 1 . This result was then confirmed by numerical simulations of 
poly tropic turbulent flows (Fig.|3] right panel). 
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Fig. 3 Left: Lognormal density PDFs for isothermal one-dimensional simulations at various Mach 
numbers, indicated by the labels. The independent variable is i = Inp. Right: Density PDFs for 
polytropic cases, with effective pol ytropic exponent y^. = 0.3 (top) and y^ = 1.7 (bottom). From 
iPassot & Vazguez-Semadenil jl998h . 



4 Turbulence and Thermodynamics 



In the previous sections we have separately discussed two different kinds of physical 
processes operating in the ISM: radiative heating and cooling (to which we refer col- 
lectively as net cooling, nA — F), and compressible turbulence. However, since both 
operate simultaneously, it is important to understand how they interact with each 
other, especially because the mean density in the Solar neighborhood, (n) ^ 1 cm^^, 
falls precisely in the thermally unst able range. This problem has been investigated 
numerically by various groups (e.g.jHennebelle & Perault),'l999';'Walder & Folinu 
2,000; Koyama 



various groups (e.g.JHennepelle & Ferault?, 'lyyyvwalder & rolini 
& Inutsuk a. 20 001 120021: IVazquez-Semadeni et al, 2000a, .2003 
2006|: [Gazol et al.[ boOli |2005l jKritsuk & Normani l2002i: ISanchez-Salcedo et al 



Piontek & Os trikeV, '2004 120051: lAudit & Hennebellel 120051: iHeitsch et al 



2002; 



2005 



Hennebelle & Audit. .2007i) . 



4.1 Density and Pressure Distributions in the ISM 

The main parameter controlling the interaction between turbulence and net cooling 
is the ratio rj = Tc/Tt, where Tc ~ e / (nmupA) is the cooling time and Tt « L/U 
is the turbulent crossing time. The remaining symbols have been defined above. 
In the limit rj ^ I, the turbulent compressions evolve much more rapidly than 
they can cool, and therefore behave nearly adiabatically. Conversely, in the limit 
T] <C 1, the fluctuations cool down essentially instantaneously while the turbu- 
lent compression is evolving, and thus they tend to reach the thermal equilibrium 
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pressure fig as soon a s they are producecj^ ( Elmegreen , 1991 : Passot et al. , 1995 



Sanchez-Salcedo et al. , 2002 : Vazquez-Semadeni et al. , 20031 : Gazol et al. , 20051) . 



Because in a turbulent flow velocity fluctuations of a wide range of amplitudes and 
size scales are present, the resulting density fluctuations in general span the whole 
range between those limits, and the actual thermal pressure of a fluid parcel is not 
uniquely determined by its density, but rather depends on the details of the velocity 
fluctuation that produced it. This causes a scatter in the values of the pressure around 
the thermal-equilibrium value in the pressure-density diagram (Fig. ID left panel), 
and also produces significant amounts of gas (up to nearly half of the total mass) 
with densities and t e mperatures in the classically forbidden thermally unstable range 
dGazol et al.L I2OOII: Ide Avillez & Breitschwerdli I2OO5I I Audit & Hennebell3. 12005 



Mac Low et al.L l2005h. a result that has been e ncountered by various observational 



studies (e.g.. lDickev et al.Lll978l : lHeilesll2001l) . In any case, the tendency of the gas 



to settle in the stable phases still shows up as a multimodality of the density PDF, 
which becomes less pronounced as the rms turbulent velocity increases (Fig.|4] right 
panel). 




Loq n (cm^^) -4-20 2 4 

^ ^ ' log n (cm-') 



Fig. 4 Left: Two-dimensional histogram of the grid cells in the pressure-density diagram for a two- 
dimensional simulation of turbulence in the thermally-bistable atomic medium, with nns velocity 
dispersion of 9 km s^' , a numerical box size of 100 pc, and the turbulent driving applied at a scale 
of 50 pc. Right: Density PDF in simulations like the one on the left panel, but with three different 
values of th e rms velocity: 4.5 km s^' (solid line), 9 km s^' (dotted line), and 1 1.3 km s^' (dashed 
line). From lGazol et all KOOSh . 



* Note that it is often believed that fast cooling implies isothermality. However, this is an en'oneous 
notion. While it is true that fast cooling is a necessary condition for isothermal behavior, the reverse 
implication does not hold. Fast cooling only implies an approach to the thermal equilibrium con- 
dition, but this need not be isothermal. The precise form of the effective equation of state depends 
on the details of the functional dependence of A and F on T and p . 
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4.2 The formation of dense, cold clouds 

Another important consequence of the interaction of "turbulence" (or, more gener- 
ally, large-scale coherent motion s of any kind) and TI i s that t he former may non- 



linearly induce the latter Indeed, iHennebelle & Peraultl(ll999[ see also Koyama & 



Inutsuka 2000) showed that transonic (i.e., with Ms ~ 1) compressions in the WNM 
can compress the medium and bring it sufficiently far from thermal equilibrium that 
it can then undergo a phase transition to the CNM (Fig.|5] left panel). This process 
amounts then to producing a cloud with a density up to lOOx larger than that of 
the WNM by means of only moderate compressions. This is in stark contrast with 
the process of producing density fluctuations by pure supersonic compressions in, 
say, an isothermal medium, in which such density contrasts would require Mach 
numbers Ms ~ 10. 



t- 7.7i (Myrfi) 'iz.o 




log n [cm-3] -'° ^ 

Fig. 5 Left: Evolutionary path (dashed line) in the P vs. p diagram of a fluid parcel initially in the 
WNM after suffering a transonic compression that nonlinearly triggers TI. The solid and dotted 
lines show the locus of Peq{p), the solid sections corresponding to linear stability and the dotted 
ones to linear instability. The solid section to the left of the dotted line corresponds to the WNM 
and the one at the right, to the CNM. The perturbed parcel evolves from left to right along the 
dashed line. From Koyama & Inutsukal (2000). Right: Projected (or column) density of a numerical 
simulation of the formation of a dense cloud formed by the convergence of two large-scale streams 
of WNM. The projection is along lines of sight perpendicular to the direction of comp r ession . The 
cloud is seen to have become turbulent and highly fragmented. From lHennebelle etalj l l2008h . 



Moreover, the cold clouds formed by this mechanism have typical sizes given by 
the size scale of the compressive wave in the transverse direction to the compression, 
thus avoiding the restriction of having the size scale of the fastest growing mode of 
TI, which is very small 0.1 pc; cf. Sec. El). The initial stages of this process 
may prod uce thin CNM sheets (|Vazauez-Semadenie t al., 2006), which are in fact 
observed ( Heiles & Troland , 20031) . However, such sheet s are quickly des tabilized, 
apparently by a combination of nonlinear thin shell (NTSI; ^Vishniac |, 1994 ), Kelvin- 
Helmholtz and Rayleigh-Taylor instabilities dHeitsch et al. . 2005), fragmenting and 
becoming turbulent. This causes the clouds to become a complex mixture of cold 
and warm gas, where the cold gas is distributed in an intrincate network of sheets. 
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filaments and clumps, possibly permeated by a dilute, warm background. An exam- 
ple of this kind of structure is shown in the right panel of Fig.|5] 



5 Turbulence in the Ionized ISM 



As discussed in the previous sections, the ionized and atomic components of the 
ISM consist of gas in a wide range of temperatures, from T ^ 1 0^ K for the HIM, to 
r ~ 40 K for the CNM. In particular. .Heiles & Troland I (l2003b report temperatures 
in the range 500 < T < 10"^ K for the WNM, and in the range 10 < T < 200 K 
for the CNM. The WIM is expected to have T ^ 10'^ K. Additionally, those same 
authors report column density-weighted nns velocity dispersions a,. 11 km s^' 
for the WNM, andgy ~7kms~' for th e CNM. Since the adiabatic sound speed is 
given by (e.g., lLandau & Lif shitzl [l959l) 



ljkT_ 



10.4 km s" 



I lo^Ky 



1/2 



(10) 



it is clear that the warm, or diffuse, gas is transonic (Ms ~ 1), while the cold, or 
dense, gas is strongly supersonic (3 ^ Ms ^ 20). Indeed, collecting measurements of 
interstellar scintillation (fluctuations in amplitude and phase of ra dio waves caused 



by sca ttering in the ionized ISM) from a variety of observations, [Armstrong et al 



(119951) estimated the power spectrum of density fluctuations in the WIM, finding that 
it is consistent with a Kolmogoro v spectrum, characteristic of incompressible tur- 
bulence (see, e.g., the rev i ews b y IVazquez-Semadeni Il999t iMac Low & KlessenL 
20041: lElmegreen & Scalo[l2004 . on scales 10^ < L < lO'^ cm. 

Mo re recently, using data from the Wisconsin Ha Mapper Observatory, Chepurnov & LazarianI 
(l2010l) liave been able to extend this result to scales ~ lO''^ cm, suggesting that the 
WIM behaves as an incompressible turbulent flow over size scales spanning more 
than 10 orders of magnitude. Although the WIM is ionized, and thus should be 
strongly coupled to the magnetic field, the turbulence then being magnetohydrody- 
namic (MHD), Kolmogorov scalin g should still apply, accord ing to the theory of 
incompressible MHD fluctuations dGoldreich & Sridhan,ll995k . The likely sources 
of kinetic energy for these tu rbulent motions are stellar e nergy sources such as su- 
pernova explosions (see, e.g.. lMac Low & KIessenll2004l) . 
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6 "Turbulence" in the Atomic and Molecular ISM 



6.1 The atomic medium 



In contrast to the relatively clear-cut situation for the ionized ISM, the turbulence 
in the neutral (atomic and molecular) gas is more complicated, and is currently un- 
der strong debate. According to the discussion in Sec. [5] the temperatures in the 
atomic gas may span a continuous range from a few tens to several thousand de- 
grees, and have velocity dispersions of Cy ^ 7-10 km s^\ suggesting that it should 
range from mildly to strongly supersonic. But because the atomic gas is thermally 
bistable (i.e., has two stable thermodynamic phases separated by an unstable one; 
Sec. |2]i, and because transonic compressions in the WNM can nonlinearly induce 
TI and thus a phase transition to the CNM (Sec. 14. Il l, the neutral atomic medium is 
expected to consist of a complex mix ture of gas spanning over two orders of mag- 



nitude in density. Early models (e.g.. iField et al.L 11969 ; lMcKee&Ostrikeit 119771) 



proposed that the phases were completely separate, but the results reported in Sec. 
14.1 [ suggest that significant amounts of gas exist as well in the unstable range, transit- 
ing between the stable phases. Numerical simulations of such systems suggest that 
the velocity disperion within the densest "clumps" is subsonic, but that the veloc- 
ity dispersion of the clumps within the diffuse substrate is supersonic with respect 
to the clumps' sound speed (although subsoni c with respect to the warmest gas; 
Kovama & Inu tsuka. 2002; Heitsch et al., '2005'). Also, note that, contrary to earlier 



ideas (e.g. lKwark, 1979; .Blitz & Shu.. 1980) . the clumps in the modern simulations 
actually form from fragmentation of large-scale clouds formed by large-scale com- 
pressive motions in the WNM, rather than the clouds forming from random collision 
and coagulation of the clumps. The complexity of this type of structure is illustrated 
in the right panel of Fig. |5] 



6.2 The Molecular Gas 



The discussion so far, involving mainly turbulence and thermodynamics, has re- 
ferred to the ionized and atomic components of the ISM. Howeve r, molecular clouds 
(MCs) have long bee n known to be str ongly self-gravitating (e.g.. iGoldreich & Kwan , 



1974: Larson , 1981). In view of this, Goldreich & Kwan 



(119741) initially proposed 



that MCs should be in a state of gravitational collapse, and that the observed mo- 
tions in MCs (as derived by the non-thermal linewid ths of molecular lines) corre- 
sponded to this collapse. However, shortly there after. IZuckerman & Palmed (119741) 
argued against this possibility by noting that, if all the molecular gas in the Galaxy 
(^moi ~ lO'^Mfg) were in free-fall, a simple estimate of the Galaxy's star formation 
rate (SFR), given by SFR ~ Mmoi/fff - 200 M0 yx-\ where Tff = yj2,n/2,2Gp is 
the free-fall time, would exceed the observed rate of '- ^ 3 Mt^ yr~ ' by roughly tw o or- 



ders of magnitude. This prompted the suggestion ( Zuckerman & Evansi 1 1 974l) that 
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the non-thermal motions in MCs corresponded instead to small-scale (in compari- 
son to the clouds' sizes) random turbulent motions. The need for these motions to 
be confined to small scales arose from the need of the turbulent (ram) pressure to 
provide an isotropic pressure that could cou nteract the clou ds' self-gravity at large, 
maintaining them in near virial equilibrium ( Larson ^ 198lh. Because tur bulence is 
known to be a dissipative phenomenon (e.g., iLandau & Lifshitzl [l959l) . research 
then focused on finding suitable sources for driving the turbulence and avoiding 
rapid dissipatio n. The main driving source was cons i dered to be energy injectio n 
from stars (e.g.. Norman & Silkl [l980l: iMcKeel [l989HMac Low & Klessenll2004 . 
and reduction of dissipation was proposed to be accomplished by having the turbu- 
lence being M HD, and consistin g mostly of Alfven waves, which do not dissipate 



as rapidly (e.g.. lShu et al.Lll987h 



However, in the last decade several results have challenged the turbulent pressure- 
support scenario: 1) Turbulence is known to be characterized by having the largest- 
velocities occur at the largest scales, and MCs are no exception, exhibiting scaling 
relations between velocity dispers ion and size suggesting that the largest dispersions 
tend t o occur at the largest scales ( Larson , 1981 ; Heyer & BruntL 2004 ; Brunt et al 



2ooi Fig. m left and middle panels). This is inconsistent with the small-scale re- 



quirement for turbulent support. 2) It was shown by several g roups that MHD turbu 



lence dissipates just as rapidly as hydrodynam ic turbulence (iMac Low et al.L 11998 



Stone et al. , 1998 ; Padoan & Nordlund 1999 ), dismissing the notion of reduced dis- 



sipation in "Alfven-wave turbulence", and thus making the presence of strong driv- 
ing sources for the turbulence an absolute necessity. 3) Clouds with very different 
contributions from various turbulence-driving mechanisms, including those with lit- 
tle or no star formation activity, such as the so-csWed Maddalena's cloud, show sim- 



ilar turbulence characteristics (IWilliams et al.L 119941 [Schneider et al.L 1201 11) . sug 



gesting that stellar energy injection is not the main source of turbulence in MCs. 

Moreover, simulations of dense cloud formation have shown that, once a large 
cold CNM cloud forms out of a collision of WNM streams, it quickly acquires a 



large enough mass t hat it can begin to co llapse gravitationally (IVazquez-Semadeni et al 



20071 l2010i I2OIII: iHeitsch et aU l2008a b). The enhancement in its column den- 



sity promotes the formation of molecular hydr ogen (H2) jHartmann et al.L 12001 



Bergin et al.L I2004L iHeitsch & HartmannL l2008h . Thus, it appears that the forma- 



tion of a molecula r cloud may require previous gravitational contraction (see 
also iMcKeei Il989l) . In addition, according to the discussion in Sees. 14.21 and 
16.11 the CNM clouds formed by converging WNM flows should be born turbu- 
lent and clumpy. The turbulent nature of the clouds furth er promotes the for- 
mation of molecular hydrogen (jOlover & Mac Low , 2007 ). The simulations by 
( Vazquez-Semadeni et al. , 2007 , 2010i.]201ll) show that the nonlinear, turbulent den- 
sity fluctuations can begin to collapse before the global collapse of the cloud is com- 
pleted (Fig. |6] right panel), both because their densities are large enough that their 
free-fall time is significantly shorter than that of the whole cloud, and because the 
free-fall time of a flattened or elongated cloud may be much larger th an that of an 



appro ximately isotropic clump within it of the same volume density (IToala et al 



20121). However, the turbulent velocities initially induced in the clouds by the con- 




Fig. 6 Left and middle panels: Second eigenimages obtained by Principal Component Analysis 
of spectroscopic data of the star-fomiing regi on NGC 1333, show ing the main contribution to the 
linewidth of molecular emission in this region jBrunt et aUl2009h . The middle image shows the re- 
gion enclosed in the rectang le in the left i mage. Black and white colors represent oppositely-signed 
components of the velocity. iBrunt et aTl j2009l) describe the pattern as a "dipole", in which large- 
scale patches of alternating velocity direction are observed. This is seen in both the large-scale 
and the small-scale images. Right panel: Image of the projected density field of a 3D numerical 
simulation with cooling, self-gravity, and magnetic fields, representing the formation of a dense 
atomic cloud by the collision of WNM streams in the direction perpendicular to the plane of the 
figure. The time shown is 20 Myr after the start of the simulation. The black dots denote "sink" 
particles, which replace local collapsing zones in the simulation. The whole cloud is also collaps- 
ing, although its collapse is not completed yet by the end of the simulation, at / = 3 1 Myr. From 
[Vazauez-Semadeni et alj i201 ih . 



verging flows in flie simulations of IVazquez-Semadeni et al.l (120071) were observed 
to be smafl compared to the velocities that dev elop due to the subseq uent gravita- 
tional contraction, while in the simulations of ( Banerjee et al. . 20091) . the clumps 
with highest internal velocity dispersions were those that had already formed col- 
lapsed objects ("sink" particles), although energy feedback from the sinks was not 
included. 

All of the above evidence suggests that the observed supersonic motions in 
molecular clouds may have a significant, perhaps dominant, component of infalling 
motions, with a (possibly subdominant) superposed random (turbulent) component 



remain ing from the initial stages of the cloud (Bate et al., 2003; Ba llesteros-Pared es et al, 
2011allbl) . and perhaps somewhat amplified by the collapse dVazquez-Semadeni et all 



19981) . In this scenario of hierarchical gravitational fragmentation, the first struc- 
tures that complete their collapse are small-scale, high-amplitude density fluctu- 
ations that are embedde d within larger-scale, smaller am plitude ones, which com- 
plete their collapse later ( Vazquez-Semadeni et al. . 20091) . The main role of the truly 
turbulent (i.e., fully random) motions is to provide the nonlinear density fluctuation 
seeds that will collapse locally after the global contraction incre ases their density 
sufficiently for them to become locally gravitationally unstable dClark & Bonnell , 
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20051) . Evide nce for such multi-scale collapse has recently begun to be observation 



[laps 

aUy detected dOalvan-Madrid et al.ll2009tlschneider et alll2010h . 



7 Summary and Conclusions 

In this contribution, we have briefly reviewed the role of two fundamental physical 
processes (net radiative cooling and self-gravity) that shape the nature of interstellar 
turbulence. We first discussed the effects of the net thermal effects that arise from 
several microphysical processes in the ISM, such as the emission of radiation from 
various ions, atoms, and molecules, which caries away thermal energy previously 
stored in the particles by collisional excitation, thus cooling the gas, and photo- 
electric production of energetic electrons off dust grains by background stellar UV 
radiation, which heats the gas. The presence of radiative heating and cooling im- 
plies in general that the gas behaves in a non-isentropic (or non-adiabatic) way, and 
in particular it may become thermally unstable in certain regimes of density and 
temperature, where small (i.e., linear) perturbations can cause runaway heating or 
cooling of the gas that only stops when the gas exits that particular regime. This in 
turn causes the gas to avoid those unstable density and temperature ranges, and tend 
to settle in the stable ones, thus tending to segregate the gas into different phases 
of different densities and/or temperatures. In classical models of the ISM, only the 
stable phases were expected to exist. 

We then discussed the interaction between trans- or supersonic turbulence, which 
produces large (i.e., nonlinear) density and velocity, and thermal instability (TI). We 
first discussed the probability density function (PDF) of the density fluctuations, 
which takes a lognormal form in isothermal regimes, and develops power-law tails 
in polytropic (i.e., of the form P oc pl'e) ones. We then noted that, since turbulence 
is an inherently mixing phenomenon, it opposes the segregating effect of thermal 
instability, causing the production of gas parcels in the classically forbidden unstable 
regimes, which may add up to nearly half the mass of the ISM, although the density 
PDF in general still exhibits some multimodality due to the gas' preference to settle 
in the stable regimes. The existence of gas in the unstable ranges has been confirmed 
by numerous observational studies. 

We next discussed the nature of the turbulence in the different ranges of den- 
sity and temperature of the gas, noting that in the diffuse ionized regions, where 
the flow is transonic (i.e., with Mach numbers Ms ~ 1), the gas appears to behave 
in an essentially incompressible way, exhibiting Kolmogorov scalings over many 
orders of magnitude in length scale. However, in the neutral atomic regions, where 
the gas is expected to be thermally unstable under the so-called isobaric criterion, 
the flow is expected to exhibit large density and temperature fluctuations, by up to 
factors ~ 100, thus being highly fragmented. Numerical simulations of this process 
suggest that the gas is transonic with respect to the warm diffuse component, but 
supersonic with respect to the cold, dense one, although those supersonic motions 
appear to correspond more to the velocity dispersion of the dense clumps within 
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the warm substrate, rather than to the internal velocity dispersion within the clumps 
themselves. 

Finally, we pointed out that large-scale compressions in the warm neutral gas, 
which may be triggered by either random turbulent motions, or by yet larger-scale 
instabilities, may nonlinearly induce the formation of large regions of dense, cold 
gas; much larger, in particular, than the most unstable scales of TI, which have 
sizes ^0.1 pc. These large clouds may easily be large enough to be gmvita- 
tionally unstable, and numerical simulations of their evolution suggest that they 
rapidly engage in gravitational contraction. The latter may in fact promote the for- 
mation of molecules, so that the clouds are likely to become molecular only after 
they begin contracting. In addition, the clouds are born internally turbulent by the 
combined effect of TI and other dynamical instabilities, and the resulting nonlin- 
ear density fluctuations ("clumps") may themselves become locally gravitationally 
unstable during the contraction of the whole large-scale cloud. Because they are 
denser, they have shorter free-fall times, and can complete their local collapses 
before the global one does, thus producing a regime of hierarchical gravitational 
fragmentation, with small-scale, short-timescale collapses occurring within larger- 
scale, longer-timescale ones. It is thus quite likely that the flow regime in the dense 
molecular clouds corresponds to a dominant multi-scale gravitational contraction, 
with smaller-amplitude random ("turbulent") motions superposed on it. Interstellar 
turbulence is seen to involve an extremely rich and complex phenomenology, even 
more so than the already-fascinating regimes of terrestrial turbulence. 
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